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Abstract 

In this paper we focus on diversity-induced resonance, which was recently found in bistable, 
excitable and other physical systems. We study the appearance of this phenomenon in a purely eco- 
nomic model of cooperating and defecting agents. Agent's contribution to a public good is seen as 
a social norm. So defecting agents face a social pressure, which decreases if free-riding becomes 
widespread. In this model, diversity among agents naturally appears because of the different sensitiv- 
ity towards the social norm. We study the evolution of cooperation as a response to the social norm 
(i) for the replicator dynamics, and (ii) for the logit dynamics by means of numerical simulations. 
Diversity-induced resonance is observed as a maximum in the response of agents to changes in the 
social norm as a function of the degree of heterogeneity in the population. We provide an analytical, 
mean-field approach for the logit dynamics and find very good agreement with the simulations. From 
a socio-economic perspective, our results show that, counter-intuitively, diversity in the individual 
sensitivity to social norms may result in a society that better follows such norms as a whole, even if 
part of the population is less prone to follow them. 



Or 



1 Introduction 

The ever-increasing interest by physicists to contribute to understanding collective phenomena in so- 
cial systems HI has mostly concentrated around highly stylized models, often directly borrowed from 
physics, using vague plausibility arguments to justify their social context 0. In this paper, we follow 
a less common route, namely to work with a model which is established in, and directly taken from, 
the social sciences. It studies the effect of social norms on the emergence of cooperation. We study its 
dynamics from the physical perspective of diversity-induced resonance, to shed new light on sustainable 
cooperation in a society where some fractions do not adhere to support it. 

In a system consisting of distinct and non-identical elements, diversity-induced resonance can be defined 
as the appearance of a maximum response to an external signal, dependent on the degree of diversity. 
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This phenomenon was first reported in Ref. [3], in the context of coupled bistable or excitable systems 
that are subject to a sub-threshold signal. It was shown that there is an optimum level of the diversity 
(quenched noise) of the coupled units that maximizes the response to the signal. Subsequent works EJE] 
[7J [3 [lOl \TT\ [121 [T31 showed that similar behavior can be observed in other physical systems, thus 
reinforcing the notion that this type of resonance can be a quite general phenomenon. In fact, diversity- 
induced resonance was also shown to appear in models related to sociophysics: It was found in discrete 
models of opinion formation lfT4l -such as the Galam model |[T5l (related to the random-field Ising 
model at zero temperature iPToll - and in continuous ones [17], of which the Deffuant model lfT8l is a 
paradigmatic example. In all cases, the average opinion synchronizes to external signals or influences 
when the diversity in the preferred opinions attains an optimum value. In a broader context, diversity- 
induced resonance can be generalized to other sources of disorder in the internal dynamics of the system 
constituents. Interestingly, even repulsive and evolving patterns of interactions can trigger a common 
collective behavior, be it synchronization |[T9l 1201 . an amplification of an external signal ifTTl I2D or 
a nonlinear increase in the volatility of the global dynamics ll22l . In a social context, these repulsive 
interactions would represent contrarians, i.e. individuals that oppose any type of consensus ll23l l24l or 
that intend to destabilize the system itself, such as the joker-like players studied in the context of social 
dilemmas E5l . 

The research reported here generalizes diversity-induced resonance by demonstrating its appearance in 
a purely economic model of social norms and their effect on cooperation |p26l . Instead of relying on a 
model rooted in physics, we study an established model from the economics literature in which diversity 
and external driving are introduced based on economic considerations. In this model, diversity appears 
naturally as an idiosyncratic propensity to follow a social norm. We demonstrate for this model that 
there is an optimal range of diversity, which leads the society to follow such norms as a whole. Different 
from the setup of diversity-induced resonance models usually studied in the physics literature, in this 
case diversity appears in a multiplicative manner and its dynamics are given by approaches typical of 
evolutionary game theory. 

The paper is organized as follows: Section |2] presents our model and its economic context. Section |3] 
summarizes our simulation results, obtained for two different types of evolutionary dynamics, to demon- 
strate the robustness of the observations. To better understand the origin of the collective dynamics, we 
present our findings for three levels of increasing modeling complexity, without and with diversity, and 
with external forcing. Subsequently, Sec. |4] improves our understanding by means of an analytical ap- 
proach for the stationary level of cooperation, whereas Sec. 14.21 investigates the response to the external 
signal. Finally, Section [5] summarizes our conclusions and discusses the implications of this work. 
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2 Model 

2.1 Economic context 

In this paper, we model conditional cooperation, a phenomenon observed in many human interactions. 
This term was introduced by Keser and van Winden E71 and Fischbacher et al. |[28l to refer to the 
fact that people often condition their cooperation on the cooperativeness of others or on their beliefs 
about others' behavior. In the specific context of Prisoner's Dilemma ||29l 1301 or Public Goods games 
11311 . this means that people are ready to contribute more to the common welfare if others contribute 
as well. Furthermore, this willingness increases with the number of contributors in the game. There is 
a large body of experimental evidence supporting the existence of this type of behavior 021 . even in 
structured populations Il33l[34ll35l . It is only consistent to ask (i) for a deeper theoretical understanding 
of these observations and (ii) their consequences for economic reasoning. The first question is partly 
answered by the theory of social preferences 11361 . that posits that non-monetary contributions to the 
utility function arise from social considerations, such as, e.g., inequity aversion or reciprocity. It has 
been argued that social preferences arise through social norms, i.e., rules of conduct that are enforced 
by internal or external sanctions [37]. Explanations for the emergence and robustness of such norms in 
evolutionary terms have been advanced ll38l l39l 1261 . thus closing the rationale to explain conditional 
cooperation in terms of social preferences. 

In this paper, we focus on the issue of norms and on the consequences of having a diverse population of 
conditional cooperators interacting in a Public Goods setup. Thus, we investigate how diversity influences 
the response to exogenous efforts to promote cooperation through social norms. Following Spichtig and 
Traxler Coll . we consider that a norm against free-riding is enforced (internally or externally). This is 
achieved by adding a contribution to the utility function such that free-riding (i.e., not contributing to 
the public good while benefiting from it) is heavily punished when rare, but the punishment weakens 
as free-riding becomes more abundant in the population. This norm leads to conditional cooperation 
because of more willingness to cooperate when the population is mostly cooperative, and the propensity 
to cooperate decreases if less participants cooperate. 

In the above context, we address the following question: How does the behavior of the population change 
if the social norm responsible for establishing a conditionally cooperative strategy varies in time? This 
question is important for two reasons. First, social norms are known to change in time, endogenously or 
exogenously, in periodic or random manners l40ll4D . Therefore, it is most important to understand how 
those changes affect the observed behavior in order to assess the stability of cooperative environments. 
Second, understanding the response of the population to changes in the current social norm can help 
policy makers to design incentives or new norms that lead to more cooperative outcomes. However, 
it should also be realized that the effort of steering the norms towards a preferred direction is costly 
and, at some point, the benefit of improving the behavior of the populations may be lower than that of 
continuing changing the norm. Therefore, assessing the optimum amount of effort invested in modifying 
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a given norm is a very relevant issue. Finally, we will come to the issue of diversity-induced resonance by 
considering that the sensitivity to the social norm depends on the individual through a specific coefficient 
to be introduced in the utility function. In the following, we will show that these issues can be addressed 
-and are related to- the phenomenon of diversity induced resonance in this system 

2.2 Model definition 

Let us now implement the ideas above in a well-defined model built on the original proposal by Spichtig 
and Traxler |[26l . We consider a population of N agents which can take one of two possible (opposite) 
actions cr, € {0, 1}, for i = 1 • • • N. We assume that "cooperative" agents take action a = 1, this way 
contributing to a public good, while "free riders" take action a = and do not contribute to the public 
good. Defining the density of cooperators as n c = N c /N , and the density of free-riders as rif = 1 — n c = 
Nf/N respectively, the utility (or payoff) function per agent is defined as 

N 



The first term in Eq. £T|) represents the cost c per agent for providing the public good, which applies only 
if agent i is cooperative, Uj = 1. The second term represents the benefit r/N per agent resulting from 
the public good. It applies regardless of the agent's action crj. Both terms describe the utility function of 
a classical public good game. The third term, new to the model, describes an additional effect resulting 
from the existence of a social norm, or social pressure, to cooperate. Free-riders with oi = face an 
(internal or external) sanction 15711 . which does not apply for cooperators with a% = 1. We assume that 
the strength of the social pressure s(n c ) depends on the density of cooperators. If n c is small, i.e. if free- 
riding is widespread, then agents deviating from cooperation may face weaker sanctions. Hence, s(n c ) 
is assumed to increase monotonously with n c , with s(l) > and lim nc _ i >o s(n c ) = 0. In the following, 
we simply choose a linear function s(n c ) = an c , with a > 0. 

Eventually, we consider that not all agents may be prone to social pressure in the same manner. To cope 
with this individual sensitivity to the social norm, we introduce a new variable 9 with realizations 9{ 
drawn from a probability distribution function g{9) with mean and standard deviation A9. Note that 
negative values of 9i imply a positive contribution to the perceived agent's utility by violating the social 
norm. This reflects the presence of contrarians/jokers E51 in the population that are willing to go against 
the system in order to benefit. Such individuals would more likely not contribute to the public good in 
presence of social pressure but, as we will see below, their presence turns out not to be an obstacle for 
the general population to conform to the social norm. 

With this utility function, the rational choice of an agent on what action to take depends on the density of 
free riders, n c and on her individual sensitivity, 9- L . Introducing c = c — r/N, it is easy to see that agents' 




+ (<Ti - 1) 9is(n c ) 



(1) 



AIM 
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decisions can be classified in three types: (i) agents will always cooperate, cjj = 1, if f?j > c/s(iV -1 ), 
(ii) agents will always free-ride, cjj = 0, if < c/s(l), and (iii) agents are conditional cooperators 
dependent on the density of free-riders in the population, i.e. they cooperate if c/s(N~ 1 ) > 9i> c/s(l). 
Note that because of lim ric _ s> o s(n c ) = 0, the criterion for the existence of cooperators is quite tight 
and often they will be absent from the population. Hence, the diversity in the individual sensitivity 9, 
precisely the standard deviation A9, will play an important role in deciding about the size of the three 
groups defined above. The final level of cooperation (as well as the influence of the social norm) will to 
a large extent be governed by the conditional cooperators. 

Finally, we will consider that the cooperation-fostering norm changes in time, which is modeled by 
assuming a time dependence of a — > a(t). This corresponds to a change of the slope of the social 
pressure function, representing periods in history where free-riding is less tolerated than in others, but 
it is always tolerated if widespread. If we further assume that agents can change their action depending 
on their expected utility, i.e. the density of cooperators n c has a dynamics defined like in the following 
section, the third term in Eq. CD representing the social pressure becomes (<jj — l)9ia(t)n c (t). Hence, 
we have a signal a(t) that changes over time because of external influences. In the present paper, for the 
sake of simplicity, and without altering the main results |[22l . we will consider a periodic change in the 
amplitude of the social norms. In absence of cooperators, the effect of this signal vanishes as well. The 
diversity in responding to the signal is given by the individual variables (cr, — i.e. only free-riders 
will face the social pressure, but they are prone to it in a heterogeneous manner. 

Studying the model in the setting of diversity-induced resonance allows us to use standard techniques for 
quantifying the response of the population to (for example) a change in the social pressure induced by a 
policy change. If the period of the signal a(t) is long enough, the results of a periodic forcing become 
equivalent to a one-time modification. Moreover, in contrast with previous studies of this phenomenon, 
the signal enters multiplicatively on the heterogeneous term. 

2.3 Evolutionary dynamics 

As mentioned above, we implement a dynamics that allow agents to change their actions dependent on the 
utility expected. For this dynamics, we use a standard evolutionary game-theoretical setup with one-shot 
games, i.e. agents have no memory about their previous action. We consider a well-mixed population, 
i.e. all agents interact together. This is dynamically equivalent to considering a mean-field version of 
the public goods game, already reflected in the sum term in Eq. (Q]). After each round of the game, 
agents collect their payoff and subsequently update their strategies according to two different dynamical 
rules, which we explain in detail below. From the various propositions for update rules in the literature 
H2ll43l . we have chosen (i) the replicator dynamics B4ll45l . which is widely used and has a well defined 
limit for N — > oo, namely the celebrated replicator equation R61, l47lk and (ii) the logit dynamics [48 ], 
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which allows for the possibility of errors or mistakes in choosing actions, and whose deterministic limit 
coincides with the best-response rule, widely used in economics Il49l . 

From a socio-economic context, both dynamic rules have a different interpretation. On the one hand, 
the replicator dynamics involves some degree of social interaction (the process is driven by imitation of 
successful strategies). On the other hand, the logit dynamics is simply based on strategic behavior. By 
choosing these quite different kinds of dynamics, we demonstrate the generality and the robustness of 
the results presented in this paper. 

Regarding the formal description, it is important to notice that the diversity for 9 introduced in our 
model no longer allows us to write down the macroscopic dynamics in terms of a single master equation. 
Instead, the system dynamics has to be split into the dynamics of groups of agents with the same value of 
9 . Let n ( 9 , a ) 58 be the number of agents with an individual sensitivity in the interval [9 — 59/2,9 + 59/2] 
choosing action a at time t (for simplicity, we also say agents are in state a at time t, i.e. "state" refers 
to "action"). Then, the rate equation for the density of cooperators with a sensitivity 9 is given by 



The transition rates to + (9) (w_ (9)) specify the overall transition into the state a = 1 (respectively, a = 0) 
for the two subpopulations with a given sensitivity 9, but different states. These transition rates depend 
on the dynamic rules chosen and are specified in the following. 

2.3.1 Replicator dynamics 

With this update rule, after every time step all agents revise their action simultaneously by selecting one 
neighbor at random, e.g. agent j, and comparing their own payoff m with their neighbor's payoff, Uj. If 
Ui > uj agent i keeps her action, whereas in the opposite case it adopts the action of the more successful 
agent j with a probability proportional to (uj — Ui). Replicator dynamics is purely imitative, meaning 
that actions not present currently in the system can not appear spontaneously. This in turn implies that 
states in which all agents defect or all contribute are absorbing states. In order to let the system leave 
those absorbing states, we have introduced noise to the dynamics: with a small probability e an agent can 
switch her action spontaneously at every time step. Subsequently, all payoffs are reset to zero and a new 
round of the game proceeds. 

For agents with an individual sensitivity 9, the overall transition rate towards the opposite state depends 
on the possible pairings with agents in the opposite state and equipped with individual sensitivity 9'. This 
yields 



n(0,l) 



ra(0,O)w+(0) -n(0,l)"-(0) 



(2) 




(3) 



(4) 



6E2 
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where g{9) is the distribution function of 0. The conditional transition rates u + (9\9'), uj^{9\9') are equal 
to the differences in payoff, if the payoff of the agent with 0' is larger, i.e. 



w_(0|0') 



u(6', 0) - u(0, +1), if u(9', 0) > u(0, +1); 
0, otherwise, 



and 

ju{9', +1) - u(0, 0), if u(Q', +1) > u(9, 0); 
0, otherwise. 



oj+(e\e') 

Using the utility function of our model, Eq. CD, these expressions become 



and 



u> + (9) = e 



, , , x \-9's(n c ) + c, if 9' < c/s(n c ): 
w_(0|0') = <^ V ' ' V h (5) 

lO, otherwise, 



ia\a>\ \-c + 0s(n c ), if > c/s{n c ); 
uj + (9\9 ) = < (6) 
[0, otherwise. 

Now, inserting Eqs. (O and © into Eqs. (O and (0]) choosing ^(0) to be a uniform distribution, we get 

w_(0) = e 

H'-K) + c) (7) 

0-A0 

+ r ( _c + ^(n c ))^ if 9 >c/s(n c ) (g) 
(^0, otherwise 

We emphasize that, in the presence of other distributions for the idiosyncratic term, the transition rates 
become more sophisticated, thus indicating the non-physical nature of this dynamics. 

2.3.2 Logit dynamics 

When considering bounded rational agents, economics literature often assumes that they do not imitate 
their neighbors, but follow a strategy or action that would yield the best payoff for them. In line with 
this assumption, one possible rule would be to change the action into cooperative (cj+) or defective (o;_) 
state with a transition rate 

W±W = 1 + exp [Tf3(u(9,l)-u(9,0))Y (9) 
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It is important to note that in this case the agent does not compare her payoff with that of another agent, 
but with the payoff she would obtain by using the opposite action. As there is no other agent involved, 
there is also no interaction term in the above equation, which makes the transition rates much simpler 
than in the previous case. This will be advantageous for an analytical approach as we will see below. 

The parameter ft in Eq. © quantifies the randomness in the process: When (3 is small, the agent is more 
likely to select another action at random, even if that action is not more successful. On the other hand, 
when (3 — > oo, the rule becomes deterministic, and the action that yields the maximum payoff is always 
chosen, as posited by Ellison R9l when introducing his (myopic) best response rule. 

3 Results 

3.1 Setup for computer simulations 

In order to present our results in a clear manner, we will deal first with the original model as introduced in 
|[26l . without considering diversity nor external forcing. This will be the baseline scenario against which 
we will subsequently illustrate the effects of diversity to proceed to our main result, namely the influence 
of an external driver and the concomitant appearance of diversity-induced resonance. 

As described in the preceding Section, the model has several parameters to specify. We start by measuring 
utilities as a function of the cost of contributing to the public good, i.e., by taking c = 1. For the 
multiplication factor we fixed r = 5 which, in a population of many agents, is too small to induce agents 
to contribute to the public good. Therefore, without the third term in Eq. (Q]) referring to the social norm, 
the only evolutionarily stable strategy is defection. For the population size, we have chosen N = 10 3 
agents (some runs were repeated with N = 10 4 for the sake of comparison, yielding the same results). 

Subsequently, we have chosen the following parameter values related to the social norm. The strength 
of the norm is given by the slope a which, in the absence of an external influence, is set as a constant 
a = 1, albeit changes of this parameter do not qualitatively modify our conclusions. Finally, for the 
sensitivity to the norm, we need to specify the parameters of the distribution g{6). In the following, we 
consider two cases: (a) There is a sensitivity to the social norm equal for all agents, which is given by the 
mean value of the distribution (homogeneous model). We will choose different values of 0. (b) The 
sensitivity to the social norm is different for all agents and randomly chosen from a uniform distribution 
in [0 — y/3A9, + y/3A9], where A9 is the standard deviation (heterogeneous model). Note that our 
choice allows for negative sensitivities with effects as described in Sect. [2] 

To monitor the evolution of the system, we have measured the time-dependent density of cooperators 
n c(t) = (V-W) X)i a i(t)- To determine the stationary level of cooperation, we compute the time-average 
number of cooperation, n c = (n c (t)) t . Subsequently, we also compute the second moment of n c (t), 
i.e. £ 2 = ((n c (t) — n c ) 2 ) t , which is the susceptibility of the system. 
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Replicator Logit 




Figure 1 : (upper row) Asymptotic fraction of cooperators n c dependent on the sensitivity to the social 
norm, 6 = 0, which is equal for all agents in the model without diversity, (lower row) Fluctuations 
£ around the expected action, which is free-riding (a = 0) for G < 2 and cooperation for > 2. 
(left column) Replicator dynamics. Different curves correspond to different values of e: triangles (0.01), 
diamonds (0.02), squares (0.05), circles (0.10). (right column) Logit dynamics. Different curves cor- 
respond to different values of /3: triangles (0.1), circles (1.0), diamonds (2.5), squares (10). The other 
parameters are described in the main text. In the upper right panel, the lines correspond to the analytical 
treatment, developed in section |4] 

3.2 Dynamics in the unforced model 
3.2.1 Model without diversity 

In the homogeneous model, the sensitivity to the social norm is equal for all agents, 6 = 0. Starting 
from an initial condition where half of the population acts as cooperators, and half as free-riders, Fig- 
ure Q] shows the asymptotic results of computer simulations for the two update dynamics introduced in 
Sect. 12.31 As it can be clearly seen in the upper panels, an increase in the parameter G -that controls 
the influence of the social norm- results in an increase in the density of cooperators. For the replica- 
tor dynamics, and for large values of randomness e, this effect becomes less visible as the width of the 
transition increases. The results for the logit dynamics point in the same direction, with /3" 1 being the 
parameter that controls the randomness or the frequency of mistakes. 

The results of both dynamics become very similar when noise is very small. Note that Fig. [T]is obtained 
for equal initial densities of contributors and free riders, but extensive simulations show that the value 
of Q at which the transition occurs does not depend on the initial condition. It is interesting to note the 
peak of the susceptibility (lower panels) close to the transition towards cooperation, both for replicator 
and logit dynamics. This is reminiscent of bistable systems which change their stability at the transition. 
It is the archetypal situation where diversity-induced resonance has already been demonstrated and, as 
we will see below, it will give rise to the same behavior in this socio-economic context. 




C.J. Tessone, A. Sanchez, F. Schweitzer: 
Diversity-induced resonance in the response to social norms 
Submitted (2012). 



Replicator Logit 




Figure 2: (upper row) Asymptotic fraction of cooperators n c and (lower row) Susceptibility £ 2 depen- 
dent on the sensitivity to the social norm, 9, different for all agents in the model with diversity. A9 
is the variance of the distribution g{6) with mean value 6 = 2. (left column) Replicator dynamics. 
Different curves correspond to different values of e: circles (0.02), squares (0.05), diamonds (0.07), 
triangles (0.10) (right column) Logit dynamics. Different curves correspond to different values of j3: 
circles (2.0), squares (2.25), diamonds (2.50), triangles (2.75). In the upper panels, we have selected 
two different initial conditions (n c (0) = 0.1 and 0.9) for both kinds of dynamics. The other parameters 
are described in the main text. In the upper right panel, the lines correspond to the analytical treatment, 
developed in section [4] 

3.2.2 Model with diversity 

Using the results from the model without diversity as a reference case, we now focus on the role of 
diversity in the sensitivity to the social norm. That means that instead of a fixed value 6 we consider an 
individual value for each agent which is drawn from the uniform distribution g{9) specified in Sect. 13. II 
The standard deviation A6> varies the degree of diversity. The results of computer simulations are shown 
in Figure From the previous discussion (cf. Fig.Q} we know that, for the chosen set of parameters, the 
transition from free riding to cooperation occurs at a value = 2. Therefore, in all the curves of Figure 
12 we have fixed the average sensitivity to this value, in order to investigate the role of diversity. When 
plotting the stationary number of cooperators, the two curves for the same parameter set correspond 
to different initial conditions with a majority of cooperators or defectors. From the simulation results, 
we can clearly conclude that diversity alone does not favor the transition towards cooperation. Also, 
increasing noise does not enhance this situation. From Fig. [2] we see that, for the replicator dynamics, the 
lower the noise the lower the cooperation in the asymptotic state, reaching the random level of n c = 0.5 
for very high values (agents make mistakes every other time step on average). For logit dynamics the 
results are similar, but for low noise we observe an asymmetric bifurcation in which the stationary state 
of low cooperation merges onto the n c = 0.5 state only for very large values of the diversity; higher 
noise values change the bifurcation toward a more symmetric form. 
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Figure 3: Response of the system in presence of a periodic square-wave forcing with logit dynamics. 
In all the plots, /3 = 2.5. Left column, first row: spectral amplification factor R, Eq. (fTTT) . Left column, 
middle row, maximum and minimum levels of cooperation attained during the evolution of n c for the 
system. Left column, lowermost row, the susceptibility. Each symbol corresponds to a different signal 
amplitude: Aa = 0.05, 0.1, 0.2, 0.5 (circle, square, diamond and triangle symbols respectively). An- 
alytical results (see main text) are represented with solid lines. In the right column, we depict the time 
dependency of the macroscopic state n c (solid, black lines), for three different values of the parameter 
AO. The values are AO = 0.7, 1.2, 1.7 in the upper, middle and lower plot, respectively. The dotted 
line represents the social pressure, while the thin line (green on-line, only in the middle plot) shows the 
signal applied (not in the same scale, for clarity). Other parameters are: T = 10 3 , N = 10 4 , r = 5, 
G = 2, a = 1. 



3.3 Dynamics under driving 

So far we have only discussed the role of the idiosyncratic sensitivity to the social norm and found that 
it does not induce a transition to cooperation. Now, as an important new ingredient, we consider that the 
influence of the norm changes in time, expressed by the time-dependent parameter a(t). Basically, any 
time dependence can be considered. For simplicity we have chosen a periodic function in the form of a 
square wave defined as 

(a + Aa, if 2nT < t < (2n + 1)T; 

at = s (1°) 
[a-Aa, if (2n + 1)T < t < 2(n + 1)T, 

with n = 0, 1, 2, .... In an adiabatic limit, where the period T is large such that the system reaches the 
stationary equilibrium in a period, this situation is equivalent to the application of a single change in the 
social pressure as perceived by individuals (by external means, like a change of policy, for example). We 
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Figure 4: Response of the system in presence of a periodic forcing with logit dynamics. In all the plots, 
Aa = 0.05. First row, spectral amplification factor R, Eq. (TTTT t. Middle row, maximum and min- 
imum levels of cooperation attained during the evolution of n c for the system. Lowermost row, the 
system susceptibility. Each symbol (color) corresponds to a different value of the inverse randomness: 
(3 = 2, 2.32, 2.5, 2.75 (circles, squares, diamonds and triangles respectively). With symbols we repre- 
sent the results obtained by means of computer simulations, while the analytical results are presented 
with solid lines. Other parameters are: T = 10 3 , N = 10 4 , r = 5, = 2, a = 1. 



have verified that using a sinusoidal function basically leads to the same results, qualitatively, than those 
shown in this Paper. So, we will focus on the expression given in Eq. (flOl ). 

We already defined the global density of cooperators n c (t) to be used as the order parameter. In particular, 
in the following, we will instead plot both the minimum and maximum values reached by the density of 
cooperators over time. To further quantify the collective response of the system to the externally changing 
influence of the social norm, we introduce the Spectral Amplification Factor (SAF), R, defined as [501 

R = A lX cV ' , — (11) 
Aa z 

Now, in addition to the variance AO of the sensitivity to the social norm, which describes an individual 
feature, we further have the change Aa in the social pressure caused by external influences. 
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Figure 5: Effect of a periodic signal applied to the social norm as a function of the population diversity 
A9 for the replicator dynamics. In all the plots, Aa = 0.05. First row, spectral amplification factor R, 
Eq. (fTTT) . Middle row, maximum and minimum levels of cooperation attained during the evolution of n c 
for the system. Lowermost row, the system's susceptibility. Each curve corresponds to a different value 
randomness: e = 0.04, 0.045, 0.05, 0.055 (circles, squares, diamonds and triangles respectively). Other 
parameters are: T = 10 3 , N = 10 4 , r = 5, 9 = 2, a = 1. 



A summary of our numerical results is presented in Fig. [3] The left column shows the spectral amplifica- 
tion factor R, the maximum and minimum values of cooperation n c and £ as a function of the standard 
deviation of the diversity, A9, for different values of the amplitude of the external driving, Aa. These 
results correspond solely to the logit dynamics; the results for the replicator dynamics are qualitatively 
similar to those presented in the plot and not shown. We further noticed that, for all choices of param- 
eters, the results are independent of the initial conditions. As it can be seen from the plots, for low A9, 
the response R is largely independent of A9. In this limit, it is possible to see the existence of super- 
threshold signal intensities Aa, which are those values exhibiting large oscillations in the limit A6 — > 0. 
For the parameters in the plot, this corresponds to Aa > 0.2. On the other hand, for smaller values of 
signal amplitude, we find that (in the limit of small heterogeneity) the system responds simply linearly 
to changes in the social norm. From a dynamical point of view, responses of the system to the external 
influence for low A9 are depicted in the right column, upper panel. 
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However, intermediate values of A8 do provide evidence for resonant behavior if the driving intensity is 
small: R shows a peak for values of A6 ~ 1, which becomes more noticeable for smaller signals. The 
oscillations of n c (t) are centered around 1/2, a value much larger than the one obtained for lower values 
of diversity. Moreover, the application of a one-time raising in the strength of the external signal may 
yield a non-linear response in terms of the growth of cooperating individuals. From a policy making point 
of view, this translates in low incentive costs being able to enforce the cooperative state throughout the 
population. When the driving amplitude Aa is much larger, the system may be able to follow the signal 
simply because the signal is super-threshold, and the same thing happens even in absence of diversity. 
Therefore, there is a true resonance phenomenon, which can be observed for low external signals, that 
elicits a strong response. In the middle panel of the right column, we show the dynamic response for 
a small applied signal, showing the large excursions in the number of cooperators when successively 
activating and de-activating the external signal. 

Finally, for very large values of diversity, no response to the external influence is observed. The amount of 
individuals with very heterogeneous responses to the external signal does not allow a significant portion 
of the population to react to the external signal, and the system's response becomes linear again. The 
latter result can be observed in a vanishing response R, small oscillation amplitudes in n c (the latter, in 
the lower panel of the right column of plots). In all the previous analyses, it is worth noticing that a peak 
in the susceptibility signals also the diversity-induced resonance in this system. 

Analyzing the role of noise for a fixed driving strength Aa, we observe another interesting feature of the 
dynamics under driving. Figs.|4]and[5]show for both the replicator and the logit dynamics the appearance 
of stochastic resonance |[50ll . That means, for an intermediate noise intensity (temperature or random- 
ness, in the proposed dynamics) the diversity-induced resonance peak is more clearly observed, whereas 
smaller or larger values of the randomness mostly suppress it. As with most stochastic phenomena, the 
resonant behavior is clearly marked also in the fluctuations of the system. The observation of stochastic 
resonance is remarkable because it shows up not in a physical, but in a socio-economic context. It in- 
dicates that some level of imperfections in the adoption of the better performing strategies may lead to 
larger responses to external stimuli. 

It is also worth mentioning that the presence of contrarians, i.e. of agents which defect even in presence 
of social pressure, can be beneficial for finding the diversity-induced resonance phenomenon. In some 
extreme cases, as shown in Fig.[4j the resonance peak may appear only in the presence of contrarians for 
P = 2.750. This finding is against our intuition that contrarians would hamper the adoption of a cooper- 
ative state in the system. It reminds on the positive influence of destructive agents on the emergence of 
cooperation in social dilemma situations as discussed in ll25l . where this phenomenon was termed 'the 
joker effect'. 



'Assuming a uniform distribution of the sensitivity 8 contrarians can be found in the population if \^3A9 > O. 
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4 Analytical approach 
4.1 Dynamics without driving 

To further understand the phenomenon of diversity-induced resonance in our model, we now develop 
an analytical approach that should be compared to the numerical simulations presented in the previous 
section. While the transition rates for the replicator dynamics, Eqs. (0 and ([8]), are too complicated for 
a tractable analytical approach, the situation is different for the logit dynamics. In this case, the density 
of agents with a given sensitivity 9 depends only on the total number of cooperators in the population, 
which is a macroscopic variable. Consequently, with the transition rates of Eq. (O and the equilibrium 
condition for the payoff function, Eq. £[]), we find for the transition rate towards the cooperative (and 
defective) states the following expression: 

W±W = l+exp{^[c-fl«- (12) 

From the above equation, we can trivially compute the density of cooperators by integrating over the 
complete population of individuals, 



which, by using the uniform distribution of the sensitivity 0, reduces to 

i r e+Ae i 



dO'- —— — — .„ . (14) 



2A9 J _ A0 l + exp{T/3[c-6>s(n c )]}' 
Expanding this equation, one readily obtains for the density of cooperators 

ln[l + exp{-/3[c-s(n c )(fl- Afl)]}] 

Hc = S{nc) + 2s{n c ) P A9 

In [1 + exp {-P[c + s(n c ){9 - Ac?)]}] 
2s(n c )pA9 



(15) 



This equation can be solved self-consistently to obtain the stationary value of n c . The corresponding 
results are shown as solid lines in Figs. \T\ and |2j right columns. We find a very good agreement between 
the numerical simulations and the prediction of our analytical approach, thus further supporting the 
validity of our results. 

In agreement with our discussion in Sect. 13.21 the system exhibits a pitchfork bifurcation. When increas- 
ing the control parameter A9, the solution n c = 1/2 changes its stability from unstable to stable, when 
the two branches (one with a majority of cooperators, the other with a majority of free-riders) collapse 
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in the center point. As observed in the simulations -and now confirmed by the analytical treatment-, the 
solutions are asymmetrical with respect to the stable point, with the lowest branch being less dependent 
on the value of the control parameter. 

This is key to understand the mechanism behind the diversity-induced resonance phenomenon in this 
socio-economic system: For intermediate values of the diversity AO, small perturbations are sufficient 
to overcome the separatrix, i.e. the unstable solution n c = 1/2 that divides the attractor basins of the 
two stable solutions. Thus, a signal which is usually too small to cause transitions between those states, 
can be sufficient to trigger such a transition near the bifurcation point. Farther from this critical point, a 
small signal only causes linear response of the system, around a stable fixed point. This fully confirms 
the discussion of the numerical results for the system with driving in the previous section. 



4.2 Relaxational dynamics with driving 

After considering the dynamics without driving in the previous section, we now turn to the dynamics 
with driving to better understand the response of the system to the external change of the norm. We note 
that the change in the density of cooperators after one state has been selected for update, is given by 

n c {t + St) = n c (t) + — ( ffi (t + St) - Ui{t)\ {a(t)}) . (16) 

where (•) represents the ensemble average, which is conditional on {a(t)}, i.e. all those states that did 
not change. Going over to small St = 1/N, we arrive at the continuous dynamics: 

^^ = (ai(t + 5t)\{a(t)})-n c (t). 

The expected value for the selected state at after update can be expressed as 

(<7i(t + St)\ {a(t)}) = Prob[a t (t + St) = 1]. (17) 

Without loss of generality, the probability that Oi{t + 5t) is +1, is given by (1 — Prob[+l — > 0]) + 
Prob[0 — > +1], which for this system is given by 

(a l (t + 5t)\{a(t)})= I ' d0'g{e'){l-uj-{e)+oj + {9)). (18) 



Restricting ourselves again to the particular case of the uniform distribution for 6 and logit dynamics, we 
have 

^ = /(n c ) = \ - n c (t) (19) 

In [cosh (ftc + /3s{n c )(A6 - 9))} 

4/3s(n c )A0 
In [cosh (/3c - /3s(n c )(A6 + 6))} 



+ 



Aj5s{n c )A6 
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If the external signal given by a(t) is slow enough, we can determine R by assuming that n c (t) reaches 
its stationary state fast compared to changes in a. Then, n c (t) = n*(a(t)). For a squared signal, the 
spectral amplification factor is simply given by 

7T (n*(a + Aa) - n*(a - Aa)) 2 
R{n c ) = . (20) 

For this forcing, the average number of cooperators reduces to n* = [n*(a + Aa) + n*(a — Aa)]/2. 
Then, the susceptibility can be computed as 

r-T/2 



r 1 

i 2 = \ dt (n*(a + Aa) - n c ) 2 + 
Jo 



T 

dt (n*(a — Aa) — n c ) 2 

I T/2 

from which we get for the susceptibility 

£ 2 = (n*(a + Aa) -n*(a- Aa)) 2 . (21) 

Figures [3] and |4]present a comparison between the analytical and numerical results. As with the previous 
comparisons, the match is very satisfactory. While our socio-economic model is quite different from a 
physics model, the dynamic observations have similar underlying mechanisms as known in physical sys- 
tems with diversity-induced resonance, which makes it possible to apply a standard analytical approach. 
For the replicator dynamics, we cannot apply the same techniques to calculate the observables. But the 
fact that we find in the simulations similarities between the logit dynamics, for which we have ana- 
lytical confirmation, and the replicator dynamics, allows us to conjecture similarities in the underlying 
mechanisms. 



5 Discussion and conclusion 

In this Paper, we have studied a socio-economic model of cooperation, to understand the effect of social 
pressure on the contribution to a public good [26]. We tried to point out analogies with the phenomenon 
of diversity-induced resonance in bistable physical systems reported in Ref. [3 ]. This was to show that 
methodological input from Physics can be beneficial for social sciences, in particular with respect to the 
vast knowledge about complex nonlinear dynamical systems. By adopting an already existing model, we 
avoided to impose a physics inspired toy model that may not have fitted the modeling paradigms of social 
sciences. 

Our analytical and numerical results demonstrate that our approach has been largely successful. Indeed, 
we found strong evidence of diversity-induced resonance, i.e., of the fact that the response of the system 
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to a weak external signal is stronger in a certain range of the parameters governing the disorder in the 
system. Importantly, such strong signals are sub-critical, meaning that these alone would not be able 
to drive a homogeneous system, whereas diversity on its own would lead to an undesired behavior (in 
our case, to a decrease in cooperation). Furthermore, we have pursued another analogy to a physical 
phenomenon, namely stochastic resonance |[50l . We found evidence that there is an optimal range of 
noise or randomness to obtain the response of the system to the external signal. 

It is most interesting to interpret the above results in terms of the original socio-economic model. In that 
context, diversity means different sensitivity to the influence of the social pressure towards behaving in a 
cooperative manner. If an external signal is emitted (e.g., changing laws or incentives by the government) 
that leads to changes of the social pressure, the population will follow these directions only if its corre- 
sponding sensitivity to such pressure is diverse, but not too little or too much. Homogeneous populations 
will simply ignore the new norms whereas very heterogeneous populations will end up behaving in some 
kind of "average" manner that does not follow the change. This is in agreement with the fact that strongly 
homogeneous groups, such as gangs or sects, are very insensitive to external influences trying to bring 
them to contribute to the general welfare. In an optimally diverse population, on the contrary, we would 
see that the most sensitive people would abide social pressure and start contributing to the common good, 
thus leading to an increment of the social pressure that pushes other agents, and so forth. 

In this context, it is important to stress that the phenomenon is robust against the kind of dynamics 
considered for the transition towards cooperation. This is particularly meaningful as the two cases studied 
in our paper, i.e. replicator and logit dynamics, correspond to two completely different approaches to 
decision making from the agent's viewpoint. While the former is based on a social, imitative, component, 
the second describes a purely strategic behavior, even a myopic one. Finally, we have observed that 
in some cases the required degree of heterogeneity for the appearance of the resonance leads to the 
existence of contrarian individuals in the population, which would benefit from going against the norm. 
This resembles the case of diversity-induced resonance arising from repulsive interactions and related 
results in social dilemmas, as mentioned in Sect. 12 

It is also worth noticing that the phenomenon of diversity-induced resonance only uses a weak signal to 
obtain the desired results. Strong signals would drive the population irrespective of its degree of diversity, 
but the external effort of the "driver" has to be much larger. This may be important for policy-making 
decisions where costly interventions in the society are not desirable because their benefit may in the 
end be smaller than the incurred cost. Of course, the requirement of diversity implies that these easily 
implemented policies may not be possible for all groups or societies, which in itself is another hint to 
policy makers about the need to estimate costs prior to specific interventions. It goes without saying that 
applications of these ideas in real life may need more complete models. For instance, one could think 
of endogenously generated norm changes, involving a feedback between actions and utility functions, or 
including the affective dimension of agents by considering their emotional response |[5T1l . On the other 
hand, applying these ideas to organizations may require a careful consideration of hierarchical effects 
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Il52l . Such improved models would lead to results that would be much more amenable to comparison 
with actual social group dynamics or even with specifically designed experiments, and thus contribute to 
our knowledge of the mechanics of social improvement. Work along these lines is in progress. 
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